Charge order in LuFe 2 04: an unlikely route to ferroelectricity 
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We present the refinement of the crystal structure of charge-ordered LuFe2 04, based on single- 
crystal x-ray diffraction data. The arrangement of the different Fe-valence states, determined with 
bond-valence-sum analysis, corresponds to a stacking of charged Fe bilayers, in contrast to the polar 
bilayers previously suggested. This arrangement is supported by an analysis of x-ray magnetic cir- 
cular dichroism spectra, which also evidences a strong charge-spin coupling. The non-polar bilayers 
are inconsistent with charge order based ferroelectricity. 
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Multiferroic oxides with strong magnetoelectric cou- 
pling are of high interest for potential information- 
technology applications but they are rare because 
the traditional mechanism of ferroelectricity is incompat- 
ible with magnetism Q ■ Consequently, much attention is 
focused on various unconventional mechanisms of ferro- 
electricity [1]. Of these, ferroelectricity originating from 
charge ordering (CO) is particularly intriguing because 
it potentially combines large electric polarizations with 
strong magnetoelectric coupling [3]. However, examples 
of oxides where this mechanism occurs are exceedingly 
rare and none is really well understood. 

LuFe 2 04 is often cited as the prototypical example of 
CO-based ferroelectricity (see, e.g., [1]). In this mate- 
rial, Fe valence order has been proposed to render the 
triangular Fe/O bilayers (see Fig. polar by mak- 
ing one of the two layers rich in Fe 2+ and the other in 
Fe 3+ H . Because of this new mechanism of ferroelec- 
tricity and also because of the high transition tempera- 
tures of CO (T co ~320 K) and magnetism (T«_~240 K) 
LuFe 2 04 is increasingly attracting attention [51425 1. That 
the Fe/O bilayers become polar upon CO has never been 
challenged. Symmetry-analysis of CO superstructure re- 
flections [2(| led to the proposal of an antiferroelectric 
stacking of the bilayer polarizations in the ground state, 
but did not cast into doubt the polar nature of the CO 
bilayers. Although these polar bilayers are generally ac- 
cepted in the LuFe204 literature [5j4l0j, a direct proof 
is lacking. An assumption-free experimental determina- 
tion of whether or not the CO in the Fe/O bilayers is 
polar would be crucial given the dependence of the pro- 
posed mechanism of ferroelectricity in LuFe2 04 on polar 
bilayers. 

In this letter, we present the first crystal structural 
refinement taking into account the superstructure due 
to CO in LuFe204, performed on single-crystal x-ray 
diffraction data. Identifying the positions of Fe 2+ and 
Fe 3+ valences in the structure with bond-valence-sum 
(BVS) analysis, an unexpected new CO pattern with 
charged Fe/O bilayers emerges (Fig.QJi). We also present 



x-ray magnetic circular dichroism (XMCD) measure- 
ments, which link (Fig, [jja,) the CO with the spin order 
determined elsewhere [23| , further corroborating the new 
CO pattern. This new CO arrangement with charged 
and non-polar bilayers is in strong contrast to all pre- 
viously suggested CO configurations with polar bilayers 
[H [2(| . We discuss the implications of this result on 
"ferroelectricity from CO" in LuFe204, addressing the 
possibility of polarizing the bilayers by an electric field. 
Finally, we also address the relevance of the strict spin- 
charge coupling to the CO transition. 

Laboratory x-ray diffraction work was done on well- 
characterized crystals with an Agilent-Technologies Su- 
perNova diffractometer using Mo-K a radiation and a 
cryojet HT for temperature control. Above Tco ( a t 
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FIG. 1. (Color online) (a) Monoclinic crystal structure C2/m 
of LuFe 2 4 (a = 5.95A, 6= 10.30 A, c= 16.96 A, /3 = 96.72°). 
The refined data was measured at 210 K. The ferrimagnetic 
high-field spin order and Fe 3+ ^ 2+ charge order is represented 
by arrows and different colors respectively, (b) Lu and O 
atoms drawn as thermal ellipsoids in projection along a. For 
comparison the Lu positions at 350 K are displayed as spheres. 
(c,d) O coordination at 210 K and 350 K for Fe 2+/3+ minority 
(red spheres indicate the O-positions at 350 K). 



2 



TABLE I. Valences from bond-valence-sum for different Fe- 
sites at 210K [C2/m] and 350K [R3m]. 
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350 K) the crystal structure of LuFe204 was refined [27 1 
in R3m symmetry, with similar results as [28| and low 
i?-factor R[F 2 > 4ct(F 2 )] = 1.87%. As already reported in 
[2fl [26l [30l ] by cooling through Tqo strong diffuse scat- 
tering along (33^) splits into sharp CO superstructure 
reflections (Fig. Oi/b), with a small incomensurability. 
Only samples showing the best magnetic behavior, cor- 
responding to those studied in , show these sharp 
superstructure peaks already at room temperature. For 
refinements the apparent small incommensuration away 
from (\\t) and (00|) type reflections was neglected, be- 
cause it most likely corresponds not to a "truly incom- 
mensurate" structure [29|, but rather to a discommen- 
suration from anti-phase-boundaries as previously pro- 
posed for LuFe204 [2(J [26| and also observed in other 
CO oxides (3l|. The superstructure reflections originate 
from three individual CO domains [20] corresponding 



to 120°-twinning with (§j|§) and symmetry-equivalent 
|i|) propagation vectors, as illustrated in 



(III) an d 
Fig. 1 of ^ 

From symmetry analysis in the hexagonal cell with 
propagation, two irreducible representations are 
allowed, both of which lower the space group to C2/m. 
These correspond to different origin positions (centers of 
inversion) of the monoclinic cell. In one case, it is located 
at the Lu positions between the bilayers, this structure 
corresponds to antiferroelectrically (AFE) stacked polar 
bilayers, as proposed in .20]. For the other case, the in- 
version center is located between the two Fe-layers of a 
bilayer, corresponding to (non-polar) bilayers with a net 
charge. This latter case was appraised as unli kely due to 
the necessity of inter-bilayer charge transfer [20( . How- 
ever, only a full structural refinement can decide which 
representation is actually realized. A consequence from 
the domain structure is that some reflections totally over- 
lap in reciprocal space, while others are difficult to sepa- 
rate, making a reliable refinement difficult. 

Therefore, a quantity of small crystals, obtained from 
one crushed sample from the same batch as in 2ol- 24j 
showing the best magnetic behavior [23}, was screened 
for their domain populations. In all experiments, the 
three domains were readily identified by the diffractome- 
ter software as twinned monoclinic cells with C2 /m sym- 
metry. Most crystals show near-equilibrium populations 
(e.g. crystal 2 in Fig.^), but some are close to a single- 
domain state (ratios of 0.03:0.09:1 for crystal 1), allevi- 
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FIG. 2. (Color online) (a) Composite precession image of 
crystal 1 in the (0 k ^Mon-plane indexed in the new mono- 
clinic cell, measured at 210 K. (b) Intensity distribution along 
(0 2£)Mon for two crystals at 210 K. 



ating the structural refinement. On this crystal we col- 
lected 8556 reflections (1285 unique); all intensities were 
corrected by numerical absorption correction using in- 
dexed crystal faces. 

A refinement in the structure model with the center 
of inversion located in the Lu-layers, corresponding to 
the representation with AFE stacked bilayers [2(| , led to 
very anisotropic displacement parameters for Lu along 
the cuex-direction. This is very unlikely for the heavy Lu 
ions. A relatively poor agreement was achieved: R[F 2 > 
4d(F 2 )]~15%. 

For refinements corresponding to the second represen- 
tation with the center of inversion between the bilayers, 
a much better R[F 2 >4cr(F 2 )] =5.96% is achieved. Addi- 
tional refinements in lower-symmetry space groups, e.g. 
Cm, allowing for both CO configurations, reproduces a 
structure very close to this second, with only marginally 
improved R- values. This makes a lower symmetry than 
C2/m very unlikely. The C2/m structural solution is 
presented in Fig. [Tk. (32|. At 210 K a Lu distortion 
along CHex with an amplitude of ^0.14 A (see Fig. [Tb) is 
clearly connected to the Fe 2+ / 3+ CO involving O shifts 
on the Fe-O-Lu path, explaining the poor refinement 
with large anisotropic displacement parameters for Lu 
on high-symmetry sites (also visible as precursor effect in 
the hexagonal solution above Tco 2£L 32]). For different 
Fe-sites strong deviations for the positions of surround- 
ing O-atoms with respect to the high-T structure are 
visible in Fig. [TJi/d, indicating a separation into two Fe- 
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FIG. 3. (Color online) (a) Fe and Lu valence states for dif- 
ferent sites from the bond-valence-sum method, (b) XMCD 
spectra across the Fe L 2 /3-edge at 120 K and 260 K (inset). 
The XAS spectra with H parallel and antiparallel to the in- 
coming beam by changed photon polarization were averaged 
by subtracting them from each other [3?| . 



valence states according to the average Fe-0 bondlengths 
(Tab. [I]) . For Fe 2+ and Fe 3+ the average Fe-0 bond length 
in a trigonal bipyramid coordination should be 2.09 A 
and 1.98 A, respectively (33l |. 

To determine the valence V from different cation sites 
a BVS analysis [34{ was performed: V = J2i exp[(rio — 
di)/0.37]. Here, di are the experimental bond lengths 
to the surrounding ions and do is a tabulated empirical 
value characteristic for the cation-anion pair [34]. The 
Lu-valences from BVS calculations are very close to 3+ 
for all temperatures and sites. The results for the Fe-sites 
are shown in Tab. Q] and illustrated by different shadings 
for different Fe sites in Fig.Q]for the CO phase. The tem- 
perature dependence of the BVS (Fig.[3k) indicates below 
Tqo & n increasing valence separation upon cooling, with 
a plateau reached below 260 K. At Tlt ~170K there is 
a magneto-structural phase transition with a small split- 
ting of structural reflections [2l|. For the data at 120 K 
an abnormal increase of the ch C x 7 calculated from the 
monoclinic lattice, is observed. However, in the refine- 
ments only very subtle changes of atom positions are 
achieved, not affecting the CO configuration (120 K in 
Fig- EJ from BVS calculations, or the C2/m symmetry. 

The Fe valence separation on majority sites tends to be 
smaller than for the minority sites. The average valence 
of all Fe-sites from BVS is ~2.4, the same as the Fe-BVS 
above Tco, suggesting non-perfect ionicity of the bonds 
to the O despite of the large valence separation. The lat- 
ter is also supported by a recent resonant x-ray diffraction 
study [5^, in which a full 2+/3+ valence separation was 
deduced from the chemical shifts of the Fe i\-edge. The 
valence separation deduced from BVS analysis is consid- 
erably larger than that for other CO Fe-oxides, except 
for Fe 2 OB0 3 

Is the new CO arrangement presented here consistent 
with the magnetic structure presented in 23]? Arc the 
magnetic structures pre-determined by the CO? X-ray 
magnetic circular dichroism (XMCD) at the Fe L-edges 
is the ideal tool to address these questions. Two previous 



XMCD studies on LuFe 2 4 were reported BBS!, but 
both were performed on samples for which no long-range 
charge and magnetic order has been demonstrated. To 
test whether the strong spin-charge coupling deduced in 
0, Gjl also applies to samples exhibiting long-range spin- 
and CO, we performed XMCD at the beamlinc 4-ID-C of 
the Advanced Photon Source (APS). We used magnetic 
fields up to 4 T ||ch C x and the incoming beam and the 
total electron yield (TEY) as x-ray absorption spectra 
(XAS). The total fluorescence yield signal (30° between 
fej and chcx) is dominated by re-absorption, but confirms 
the bulk nature of our XMCD f37|. The XMCD signal 
was then calculated from the difference between the XAS 
(from TEY) for positive and negative circular polariza- 
tion and with no non-magnetic XMCD contri- 
butions [37[ . To see if there is any change in the CO con- 
figuration or structure between the two magnetic phases, 
we have done additional high-resolution x-ray diffraction 
at the beamline 6-ID-D (APS) above T LT . The diffrac- 
tion data in i?||cHex up to 2.5 T, show neither a change 
in the CO configuration nor a structural transition. 

In the high-field ferrimagnetic phase [23|, the shape 
of the XMCD spectra A/j, (Fig. [3b) is similar to the 



ones shown in [7|, ll3j. With the sum rule |38|, |40j we 
could extract from S(A/x) a similar orbital-to-spin mo- 
ment ratio of ~0.3 corresponding to an orbital magnetic 



moment of ^0.7/j,B/f.u., as previously reported [38 . 
This observed unquenched orbital moment excludes the 
possibility of Fe 2+ orbital order for the ferrimagnetic 
phase; orbital order would imply a lifting of the two- 
fold-degeneracy of the lowest crystal-field doublet, which 
is occupied by a minority spin. However, this degeneracy 
is necessary for an orbital magnetic moment [39] ] . For the 
antiferro and paramagnetic phase it can also be excluded, 
indirectly, because there is no structural component in 
the transitions involved. Due to the structural transition 
at Tlt 2l| we can not exclude long-range orbital order 



in this low temperature phase, which could be consistent 
with the observed lattice parameter changes. However, 
a detailed discussion is beyond the scope of the present 
study. 

Two prominent peaks in the L3 region of the XAS 
are readily identified as the chemically shifted Fe 2+ and 
Fe 3+ white lines HEDl. In the XMCD spectra the large 
downward-peak at the Fe 2+ position and smaller upward- 
peak at the Fe 3+ position, directly imply that the net 
moment of Fe 2+ is in field-direction and a smaller net 
moment of the Fe 3+ sites points opposite to the field. For 
the local Fe 2+ / 3+ spin configurations the model of 0,(3]; 
extracted from a similar XMCD shape, together with the 
here presented CO is consistent with the ferrimagnetic 
spinstructure of [23| . 

More important, the above implications of the XMCD 
signal, combined with the ferrimagnetic model [23! , verify 
the novel CO configuration. Given the absence of partial 
disorder only three valence specific local spin configura- 



4 



tions are possible, of which only one is consistent with 
the overall magnetic saturation moment of ^3/ZB/f-U. 



37]: All Fe 2+ as well as i of the Fe 3+ spins are aligned 



in i^-direction, | of the Fe 3+ spi ns p oint opposite to H, 
the same model as proposed in [7j,[l3j. Combining this lo- 
cal spin-charge configuration with the ferrimagnetic spin 
order |23[ directly excludes any (anti)ferroelectric model 
preserving mirror symmetry. Ignoring mirror symmetry, 
28 configurations are possible |41| , of which however only 
the one presented in Fig.QJi fits to the right intensity dis- 
tribution along the (0 2£)m O ii line (Fig. [5p); this is also 
the only one of the 28 preserving mirror symmetry. 

Furthermore, as discussed in 23] the refinement of 
spin structures can be improved by introducing differ- 
ent magnetic moments for Fe 2+ and Fe 3+ according to 
the charged-bilayer CO model, which is not the case for 
any CO with polar bilayers. This also supports the above 
analysis, though by itself the weight of this evidence is 
reduced by a similar improvement of the refinement re- 
garding a possible magnetic contamination. 

Thus, structure refinement, XMCD, and magnetic con- 
trast, all clearly identify the CO configuration with 
charged bilayers (Fig.QJi) as the correct one. This charge 
pattern is very surprising, because it requires inter- 
bilayer charge transfer. For this reason it was considered 
before only in [20j j , where it was mentioned as symmetry- 
allowed but excluded as physically unlikely. Understand- 
ing the origin of this long range (~ 6 A) charge transfer 
calls for further theoretical work. 

Importantly, this new CO structure does not contain 
polar bilayers. in contrast to what was previously pro- 
posed (e.g. @, Hoj]), casting doubt on the "ferroelectricity 
from charge ordering" scenario. How general is our re- 
sult, given the significant reported (see, e.g., 12[) samplc- 
to-sample variations? Clearly, the structure refinement 
can be expected to be representative for all samples where 
(A, i, halfinteger)-reflections are observed as main CO 



der parameter (e.g. 0, HI, El, 13" 26 L 30|), the similarity of 
observed XMCD spectra with |Vl. Il3j even suggests that 
the same basic CO configuration also applies to samples 
without long-range CO (e.g. [14]). In particular, our re- 
finement should be valid for the samples on which pyro- 
electric current measurements have been reported [J]. 

To explain the pyroelectric current measurements, 
some of us proposed that a ferroelectric CO might 
be stabilized by an electric field, though such a scenario 
seems less likely when charged bilayers have to be polar- 
ized. Indeed, a CO remaining completely robust in elec- 
tric fields has been reported by Wen et al. [HI , based on 
neutron diffraction. We have confirmed this as also valid 
for our samples by additional x-ray diffraction measure- 
ments in electric fields up to 15kVcm _1 at APS 6-ID-D 
and 6-ID-B, and thus conclude that a ferroelectric CO 
cannot be stabilized by electric fields. 

The relatively low resistivity around Tco IS" HI could 
provide an alternative explanation for the pyroelectric 



current measurements of [J|, because in the presence of 
residual conductivity non-ferroelectrics can exhibit cur- 
rents strongly resembling ferroelectric depolarization cur- 
rents, due to space-charge effects [Hf . The also observed 
3) 01 giant dielectric constants could be attributed to 
interface effects [l9[. All reported macroscopic indica- 
tions of ferroelectric behavior in LuFe2 04 are therefore 
most likely due to extrinsic effects. 

Returning to XMCD, the analysis not only shows the 
consistency of the new CO and the spin order, but also 
implies a strict coupling of these orders. Interestingly, 
XMCD spectra taken above T N in ±4 T (Fig. [3p inset) 
have a small amplitude, but indicate the same Fe 2+ / 3+ 
spin configuration as in the ferrimagnetic phase. This is 
consistent with the conclusion for H = of randomly- 
stacked bilayers that are still individually magnetically 
ordered based on diffuse magnetic scattering [23|: par- 
tial polarization by a magnetic field is then expected to 
lead to the same relative net moments on Fe 2+ and Fe 3+ , 
provided the spin-charge coupling remains. This signi- 
fies still ordered Fe-bilayers in the paramagnetic phase 
with strictly coupled charge- and spin order persisting 
well above Tjv, from susceptibility also likely above Tco- 
This suggests at high T short-range precursors with al- 
ready coupled local spin- and CO. This coupling already 
above Tco is most likely the origin of the magnetic-field 
control of charge structures reported in [§] . 

In conclusion, crystal structure refinements of charge 
ordered LuFe204 show that the Fc/O bilayers are 
charged rather than polar. This is further supported 
by an analysis of XMCD data, which also indicates a 
strict spin-charge coupling extending to the fluctuations- 
regime above the ordering temperature. The non-polar 
CO, which is not affected by electric fields, precludes CO- 
based ferroelectricity in LuFe2 04. Hence, a clear exam- 
ple of an oxide material with ferroelectricity originating 
from CO has yet to be identified. 
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